Process control using process data and yield data

ABSTRACT

A method for monitoring a manufacturing process features acquiring metrology data for semiconductor wafers at the conclusion of a final process step for the manufacturing process (“Step a”). Data is acquired for a plurality of process variables for a first process step for manufacturing semiconductor wafers (“Step b”). A first mathematical model of the first process step is created based on the metrology data and the acquired data for the plurality of process variables for the first process step (“Step c”). Steps b and c are repeated for at least a second process step for manufacturing the semiconductor wafers (“Step d”). An nth mathematical model is created based on the metrology data and the data for the plurality of process variables for each of the n process steps (&#39;Step e“). A top level mathematical model is created based on the metrology data and the models created by steps c, d and e (”Step f&#39;). The top level mathematical model of Step f is based on those process variables that have a substantial effect on the metrology data.

CROSS REFERENCE TO RELATED APPLICATIONS

The present application is a continuation of prior co-pending U.S. application Ser. No. 12/044,346 filed on Mar. 7, 2008, the entire disclosure of which is incorporated herein by reference in its entirety.

FIELD OF THE INVENTION

The invention generally relates to methods for monitoring a manufacturing process. More particularly, the invention relates to using metrology data and process data to monitor a manufacturing process.

BACKGROUND OF THE INVENTION

Historically, semiconductor device manufacturers have managed the transition to tighter process/materials specifications by depending on process tool manufacturers to design better and faster process/hardware configurations. As device geometries shrink to the nanometer scale, however, the increasing complexity of manufacturing processes has changed the landscape that must be negotiated to meet and maintain process/materials specifications.

Stand-alone control of process tools (based on equipment state data) will not maintain viable yields at 65 and 45 nm. Advanced device processing requires tool-level control based on the combination of advanced equipment control (AEC) and sensor based process control. Furthermore, tool level control alone cannot meet all of the control needs of advanced device fabrication. System-wide implementation of advanced process control (APC) that integrates AEC with e-diagnostics, Fault Detection and Classification (FDC) and predictive mathematical models of the device manufacturing process are required.

Economic as well as technological drivers exists for the move to AEC/APC. The cost of purchasing (and developing) stand-alone process tools to meet advanced production specifications is expected to be staggering. A cost effective alternative to the purchase of a new generation of process equipment exists through the combined use of add-on sensors with AEC/APC in existing (legacy) equipment. Sensor-based AEC/APC in legacy equipment can drive these tools to the tighter specifications needed for nanometer-scale device fabrication. Additional cost benefits can be realized from reductions in scrap/rework (especially in 300 mm wafer processing) and in lower test wafer use since these can be reduced or even eliminated in systems using wafer- and/or process-state data for process control. Sensor-based AEC/APC can also reduce preventive maintenance downtimes and time to process qualifications, while increasing process capabilities, device performance, yields and fabrication throughput.

A need therefore exists for improved systems and methods for monitoring manufacturing processes. A need also exists for detecting and classifying faults associated with manufacturing processes and outputs of the manufacturing processes.

SUMMARY OF THE INVENTION

In one aspect, the invention features a method for monitoring a manufacturing process. The method involves acquiring metrology data for semiconductor wafers at the conclusion of a final process step for the manufacturing process (“Step a”). The method also involves acquiring data for a plurality of process variables for a first process step for manufacturing semiconductor wafers (“Step b”). The method also involves creating a first mathematical model of the first process step based on the metrology data and the acquired data for the plurality of process variables for the first process step (“Step c”). The method also involves repeating steps b and step c for at least a second process step for manufacturing the semiconductor wafers (“Step d”). The method also involves creating an n^(th) mathematical model based on the metrology data and the data for the plurality of process variables for each of the n process steps (“Step e”). The method also involves creating a top level mathematical model based on the metrology data and the models created by steps c, d and e (“Step f”). The top level mathematical model is based on those process variables that have a substantial effect on the metrology data.

In some embodiments, a method for monitoring a manufacturing process involves comparing data for the first process step for manufacturing a semiconductor wafer to the top level model of step f to determine whether a fault occurred during the first process step. In some embodiments, the method involves comparing data for the second process step for manufacturing a semiconductor wafer to the top level model of step f to determine whether a fault occurred during the second process step. In some embodiments, the method involves comparing data for each process step for manufacturing a semiconductor wafer to the top level model of step f to determine whether a fault occurred during any of the n process steps.

In some embodiments, the method for monitoring a manufacturing process involves comparing data for the first process step for manufacturing a semiconductor wafer to the top level model of step f to determine whether manufacturing yield will be above a predetermined limit.

In some embodiments, the mathematical model of step c is created using a multivariate analysis technique. The mathematical models of the n process steps can be partial least squares models. In some embodiments, creating the top level model involves combining scores for the partial least squares models of the n process steps and the mathematical model of the metrology data.

In some embodiments, the method involves comparing data for at least one of the n process steps for subsequent manufacturing runs to the top level model of step f to generate a prediction of yield for the manufacturing process. In some embodiments, the metrology data includes yield measurements from semiconductor wafers produced by one or more manufacturing runs.

In some embodiments, the method involves comparing data for at least one of the n process steps for subsequent manufacturing runs to the top level model of step f to generate a prediction of the metrology results for the manufacturing process.

In some embodiments, the method of monitoring a manufacturing process also involves modifying operating properties of a process step based on whether the top level model determines that a fault has occurred during the process step. In some embodiments, the mathematical model of at least one of the n process steps is updated by adjusting a weighting in the mathematical model associated with the plurality of process variables based on which variables are determined to contribute more significantly to a fault occurring during at least one of the n process steps. In some embodiments, the process variables that have a substantial effect on the metrology data are weighted more in creating the top level mathematical model than process variables that do not have a substantial effect on the metrology data. In some embodiments, the models created by steps c., d. and e. comprise process information and noise information. In some embodiments, the top level mathematical model comprises process variables that affect yield of the manufacturing process.

In another aspect, the invention features a system for creating a mathematical model according to a multivariate analysis. The system includes means for acquiring metrology data for semiconductor wafers at the conclusion of a final process step for a manufacturing process and means for acquiring data for a plurality of process variables for a first process step for manufacturing semiconductor wafers. The system also includes a model generation module for creating a first mathematical model of the first process step based on the metrology data and the acquired data for the plurality of process variables for the first process step (“Element c”). The system also includes means for acquiring data for a plurality of process variables for at least a second process step for manufacturing the semiconductor wafers and a model generation module for creating a second mathematical model of the first process step based on the metrology data and the acquired data for the plurality of process variables for the first process step (“Element d”). The system also includes a model generation module for creating an n^(th) mathematical model based on the metrology data and the data for the plurality of process variables for each of the n process steps (“Element e”). The system also includes a model generation module for creating a top level mathematical model based on the metrology data and the models created by elements c, d and e (“Element f”). The top level mathematical model is based on those process variables that have a substantial effect on the metrology data.

In some embodiments, the system for creating a mathematical model according to a multivariate analysis includes a processor for comparing data for each process step for manufacturing a semiconductor wafer to the top level model generated by element f. to determine whether a fault occurred during any of the n process steps.

In some embodiments, the system includes a processor output connected to at least one tool used to perform at least one process step. The processor can modify an operating property of at least one process step based on whether the top level model determines that a fault has occurred during at least one process step. The processor can identify which of the process variables causes manufacturing yield to be below a predetermined limit.

In some embodiments, the system includes a processor for comparing data for the first process step for manufacturing a semiconductor wafer to the top level model generated by element f to determine whether manufacturing yield will be above a predetermined limit. In some embodiments, the process variables that have a substantial effect on the metrology data are weighted more in creating the top level mathematical model than process variables that do not have a substantial effect on the metrology data. In some embodiments, the top level mathematical model is a yield prediction model based on the process variables selected by the models created by elements c., d. and e.

BRIEF DESCRIPTION OF THE DRAWINGS

The advantages of the invention described above, together with further advantages, may be better understood by referring to the following description taken in conjunction with the accompanying drawings. The drawings are not necessarily to scale, emphasis instead generally being placed upon illustrating the principles of the invention.

FIG. 1 is a schematic drawing depicting a model for monitoring a manufacturing process, according to an illustrative embodiment of the invention.

FIG. 2 is a schematic drawing depicting implementing a model for monitoring a manufacturing process, according to an illustrative embodiment of the invention.

FIG. 3 is a graph comparing a predicted yield of a manufacturing process based on a model to an observed yield of the manufacturing process, according to an illustrative embodiment of the invention.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 1 is a schematic drawing depicting a method 100 for creating a model for monitoring a manufacturing process for manufacturing a lot of semiconductor wafers 104, according to an illustrative embodiment of the invention. The method 100 involves combining process data and yield data (e.g., metrology data) to create a model for process control of a manufacturing process and yield prediction for a manufacturing process. The manufacturing process includes a plurality of process steps 108A, 108B, 108C or 108D for processing the lot of semiconductor wafers 104. In this embodiment, the manufacturing process includes “n” number of process steps 108A, 108B, 108C or 108D. The method 100 involves acquiring metrology data 112 taken from the final output 134 of the plurality of process steps 108A, 108B, 108C or 108D. In some embodiments, the step 112 involves acquiring data for a plurality of process variables for the output 134. The data acquired for each of the process steps 108A, 108B, 108C or 108D can include, for example, data from a series of semiconductor deposition process steps, ion implant steps, material removal steps performed to manufacture the final semiconductor wafers 130.

The method 100 creates models based on the plurality of process steps 108A, 108B, 108C and 108D. The method 100 involves creating a first mathematical model 116A (generally 116) from data acquired from process step 108A and the metrology data 112 acquired from the final semiconductor wafers 130. This model generation (step 116) is repeated for at least a second process step (108B). In this embodiment, generation step 116 is repeated for each process step 108, thereby generating “n” mathematical models. In this embodiment, the step 116 of creating mathematical models (116A, 116B, 116C or 116D) includes using data acquired from the plurality of process steps 108A, 108B, 108C or 108D and the metrology data 112 to generate batch level Partial Least Squares (PLS) models for each process step. The method 100 also includes the step 120 of exporting batch level scores from each model 116A, 116B, 116C or 116D and merging the batch level scores with the metrology data 112. The method 100 also includes generating a top level mathematical model 124 that is based on those process variables that have a substantial effect on the metrology data 112.

In some embodiments, creating the mathematical models based on the data acquired involves using multivariate analysis techniques or neural network modeling techniques. By way of example, multivariate analysis techniques may be implemented by using the SIMCA-P+ multivariate analysis software package (Umetrics, offices in Kinnelon, N.J.). In one embodiment, the model creation step (e.g., step 116A) involves formatting in a matrix X (see EQN. 1) the metrology data acquired, step 112 for each process step, where the matrix X has N rows and K columns. K is the number of manufacturing related variables and N is the number of outputs (e.g., observations, number of wafers) of the specific manufacturing process (e.g., process step 108A).

$\begin{matrix} {X = \begin{bmatrix} 1 & \ldots & \ldots & \ldots & K \\ \ldots & \ldots & \ldots & \ldots & \ldots \\ \ldots & \ldots & \ldots & \ldots & \ldots \\ \ldots & \ldots & \ldots & \ldots & \ldots \\ N & \ldots & \ldots & \ldots & \ldots \end{bmatrix}} & {{EQN}.\mspace{14mu} 1} \end{matrix}$

The model creation step involves calculating the average and standard deviation for each column K of the matrix X. The model creation step then involves centering and scaling each column K of the matrix X. The average of each column is subtracted from each entry in the corresponding column of the matrix X (centering). Each entry in the matrix X is then divided by the standard deviation of the corresponding column of the matrix X (scaling).

PCA Method

The Principal Component Analysis (PCA) method is then employed in developing a model. The PCA method reduces the number of dimensions (K variables) of the original data (original matrix X) down to a few independent variables that are a linear combination of the original data (K variables of original matrix X).

In one embodiment, the model developing step then involves determining the principal components and loadings of the model using an non-linear iterative non-linear partial least squares (NIPALS) method; as given by the following equation:

X=TP′+E   EQN 2

where T is the scores matrix, P is the loading matrix, P′ is the transpose of the P matrix, and E is the residual matrix. The scores matrix T is given by:

T=[t₁ t₂ . . . t_(A)]  EQN. 3

and the loading matrix P is given by:

P=[p₁ p₂ . . . p_(A)]  EQN. 4

where, for example, t₁ is the 1^(st) principal component and p₁ is the 1^(st) load vector. The 1^(st) principal component (t₁) is also referred to as the 1^(st) score vector of the X matrix and is a linear combination of the manufacturing related variables. The 2^(nd) principal component (t₂) is also referred to as the 2^(nd) score vector of the X matrix and is a linear combination of the manufacturing related variables. The subscript A is the total number of principal components used to produce the model.

The NIPALS method then involves iteratively employing the following equations:

$\begin{matrix} {p_{1}^{\prime} = \frac{t_{1}^{\prime}X}{t_{1}^{\prime}t_{1}}} & {{EQN}.\mspace{14mu} 5} \end{matrix}$

where starting with an initial condition for the t₁ vector initially equal to the column of the X matrix (EQN. 1) having the maximum variance.

The NIPALS method then involves normalizing p₁ as follows:

∥p₁∥=1.0   EQN. 6

and then employing the following equation:

$\begin{matrix} {t_{1} = \frac{{Xp}_{1}}{p_{1}^{\prime}p_{1}}} & {{EQN}.\mspace{14mu} 7} \end{matrix}$

The calculations of EQNS. 5-7 are then repeated until the calculations converge. By way of example, a user may specify (based on user experience) that the calculations have converged when

$\begin{matrix} {\frac{{t_{1} - t_{1{old}}}}{t_{1}}{\langle{{1e} - 6}}} & {{EQN}.\mspace{14mu} 8} \end{matrix}$

where t_(1old) is the value of t₁ from the prior iteration of performing the calculations using EQNS. 5-7.

The residual matrix E is then determined from the following relationship:

E=X−t ₁ p ₁ ^(t)   EQN. 9

The residual matrix E is then used in place of matrix X in EQNS. 5 and 7 to calculate the second principal component (t₂) and the second load vector (p₂). The steps associated with EQNS. 5-9 are then similarly performed for each additional principal component (t_(A)). The methods described herein are applicable to applications involving multiple numbers of principal components.

Principal component coefficients for each observation (e.g., output of the manufacturing process) are then calculated using the equations (t_(A) vectors) of the principal components (EQN. 3). Users generally use the principal component coefficients to determine if there is any clustering associated with the coefficients. Those principal component coefficients that are close to each other when plotted in A-dimensional space may sometimes have similar properties (where A is the total number of principal components used to produce the model). However, in many cases it is not possible to identify similarities and dissimilarities based on the principal component coefficients because all the original manufacturing related variables have been used to develop the model and each of the variables has been given similar weight (the principal component score is the linear combination of all the variables multiplied by its weight). Accordingly, it is therefore necessary to identify those manufacturing related variables that have a substantial effect on the model for detecting faulty outputs of the manufacturing process.

Identifying the manufacturing related variables that have a substantial effect can involve performing a partial least squares discrimination analysis (PLSDA) in which those manufacturing related variables are selected that have a substantial effect on the model for detecting faulty outputs of the manufacturing process.

PLS and PLSDA Method

The PLS and PLSDA method is employed as given by the following equations:

X=TP′+E   EQN. 10

Y=UC′+F   EQN. 11

T=XW*   EQN. 12

Y=BX+E   EQN. 13

B=W(P′W)⁻¹ C′  EQN. 14

where Y is an (M×N) matrix that contains data associated with M features of the N observations, B is the regression coefficient, and where

W*=W(P′W)⁻¹   EQN. 15

The M features may, for example, be a specific metrology metric value associated with each observation (e.g., each wafer) or a class identifier.

T is the score matrix of X, and U is the score matrix of Y. W is the X-weight matrix and C is the Y-weight matrix. F is the residual matrix associated with Y. The W matrix is given by:

W=[w₁ w₂ . . . w_(A)]  EQN. 16

and the C matrix is given by:

C=[c₁ c₂ . . . c_(A)]  EQN. 17

Similarly as previously described herein, the NIPALS method is used in accordance with the following equations:

$\begin{matrix} {w_{1} = \frac{X^{T}u_{1}}{u_{1}^{T}u_{1}}} & {{EQN}.\mspace{14mu} 18} \end{matrix}$

Where u₁ is set to be equal to one column of the Y matrix as an initial condition. The vector w₁ is then scaled to unit length by:

∥w₁∥=1.0   EQN. 19

The next step in the NIPALS method involves:

$\begin{matrix} {t_{1} = {Xw}_{1}} & {{EQN}.\mspace{14mu} 20} \\ {c_{1} = \frac{Y^{T}t_{1}}{t_{1}^{T}t_{1}}} & {{EQN}.\mspace{14mu} 21} \\ {u_{1} = \frac{{Yc}_{1}}{c_{1}^{T}c_{1}}} & {{EQN}.\mspace{14mu} 22} \end{matrix}$

where the calculations of EQNS. 18-22 are repeated until the calculations converge. By way of example, a user may specify (based on user experience) that the calculations have converged when

$\begin{matrix} {\frac{{t_{1} - t_{1{old}}}}{t_{1}}{\langle{{1e} - 6}}} & {{EQN}.\mspace{14mu} 23} \end{matrix}$

where t_(1old) is the value of t₁ from the prior iteration of performing the calculations using EQNS. 18-22.

Once the calculations have converged, the X loadings are determined from:

$\begin{matrix} {p_{1} = \frac{X^{T}t_{1}}{t_{1}^{T}t_{1}}} & {{EQN}.\mspace{14mu} 24} \end{matrix}$

The residual matrix E is then determined from the following relationship:

E=X−t ₁ p ₁ ^(t)   EQN. 25

The residual matrix F is then determined from the following relationship:

F=Y−t ₁ c ₁ ^(t)   EQN. 26

The residual matrices E and F are then used in place of matrices X and Y in EQNS. 18-26 to calculate the second principal component (t₂), second load vector (p₂), second score vector of Y (u₂), second X-weight vector (w₂), and second Y-weight vector (c₂). The steps associated with EQNS. 18-26 are then similarly performed for each additional principal component (t_(A)). The methods described herein are applicable to applications involving multiple numbers of principal components.

The next step in determining which of the K manufacturing related variables have a substantial effect on the model detecting the fault condition of the manufacturing process involves plotting the W* matrix and the C matrix. If there are two features (M), then the C matrix will have two values C₁ and C₂ for each feature which are plotted as coordinates in a plane. The W* matrix contains two values for each manufacturing related variable which also are plotted as coordinates in the same plane.

In this manner, the coordinate values of W* and C are superimposed on each other. Those manufacturing related variable coordinates (W*) that are located in the plane in proximity to a feature have a substantial effect on the model being able to predict whether a given wafer has that particular feature. Only those columns (K) corresponding to the substantial manufacturing related variables are then used to update the model using EQNS. 1-9. The updated model is then capable of distinguishing between outputs of the manufacturing process to identify which outputs are faulty.

EQNS. 1-26 can be used at an individual process level (e.g., including one or more process steps) and/or at a combined process level (e.g., including one or more processes, each of which can include one or more process steps). For example, EQNS. 1-26 as described above can be used in a method for creating a mathematical model for each process step, according to an illustrative embodiment. In one embodiment, EQNS. 1-26 as described above are used in a method for creating a mathematic model for an individual process.

Generating A Top Level Mathematical Model:

In some embodiments, an unfolding process is used to generate the top level mathematical model. Two methods can be used to generate a top level mathematical model: i) one method uses the actual data itself (e.g., data acquired for each of the process steps 108A, 108B, 108C or 108D of FIG. 1) and ii) the other method uses the scores from a lower level model. For example, EQNS. 27-29 are examples of matrix representations of the unfolded data. For example, a 3 step process can last 30 seconds with 10 wafers. In this example, the matrices for the data collection (lower level) would have the following form:

$\begin{matrix} {{{Step}\mspace{14mu} 1}:\begin{pmatrix} k_{1,1} & k_{2,1} & k_{3,1} & k_{4,1} \\ k_{1,2} & k_{2,2} & k_{3,2} & k_{4,2} \\ k_{1,3} & k_{2,3} & k_{3,3} & k_{4,3} \\ k_{1,4} & k_{2,4} & k_{3,4} & k_{4,4} \\ k_{1,5} & k_{2,5} & k_{3,5} & k_{4,5} \\ k_{1,6} & k_{2,6} & k_{3,6} & k_{4,6} \\ k_{1,7} & k_{2,7} & k_{3,7} & k_{4,7} \\ k_{1,8} & k_{2,8} & k_{3,8} & k_{4,8} \\ k_{1,9} & k_{2,9} & k_{3,9} & k_{4,9} \\ k_{1,10} & k_{2,10} & k_{3,10} & k_{4,10} \end{pmatrix}} & {{EQN}.\mspace{14mu} 27} \\ {{{Step}\mspace{14mu} 2}:\begin{pmatrix} k_{1,1} & k_{2,1} & k_{3,1} & k_{4,1} \\ k_{1,2} & k_{2,2} & k_{3,2} & k_{4,2} \\ k_{1,3} & k_{2,3} & k_{3,3} & k_{4,3} \\ k_{1,4} & k_{2,4} & k_{3,4} & k_{4,4} \\ k_{1,5} & k_{2,5} & k_{3,5} & k_{4,5} \\ k_{1,6} & k_{2,6} & k_{3,6} & k_{4,6} \\ k_{1,7} & k_{2,7} & k_{3,7} & k_{4,7} \\ k_{1,8} & k_{2,8} & k_{3,8} & k_{4,8} \\ k_{1,9} & k_{2,9} & k_{3,9} & k_{4,9} \\ k_{1,10} & k_{2,10} & k_{3,10} & k_{4,10} \end{pmatrix}} & {{EQN}.\mspace{14mu} 28} \\ {{{Step}\mspace{14mu} 3}:\begin{pmatrix} k_{1,1} & k_{2,1} & k_{3,1} & k_{4,1} \\ k_{1,2} & k_{2,2} & k_{3,2} & k_{4,2} \\ k_{1,3} & k_{2,3} & k_{3,3} & k_{4,3} \\ k_{1,4} & k_{2,4} & k_{3,4} & k_{4,4} \\ k_{1,5} & k_{2,5} & k_{3,5} & k_{4,5} \\ k_{1,6} & k_{2,6} & k_{3,6} & k_{4,6} \\ k_{1,7} & k_{2,7} & k_{3,7} & k_{4,7} \\ k_{1,8} & k_{2,8} & k_{3,8} & k_{4,8} \\ k_{1,9} & k_{2,9} & k_{3,9} & k_{4,9} \\ k_{1,10} & k_{2,10} & k_{3,10} & k_{4,10} \end{pmatrix}} & {{EQN}.\mspace{14mu} 29} \end{matrix}$

In some embodiments, the same parameters (k_(#,#)) are not measured in each process step. For example, some, all or none of the parameters can be the same for each process step. In some embodiments, each step does not have to have the same number of samples (e.g., step one could last 10 seconds, step 2 could last 20 seconds, and step three could last 2 seconds). In some embodiments, the steps do not have to be sampled at the same sampling rate. In this embodiment, EQNS. 27-29 are selected to demonstrate the unfolding process, where the form of the matrices repeats itself for each of the wafers (10 in this example) monitored.

Method 1—Unfolding of Actual Data

In this embodiment, the method involves transforming each matrix (e.g., matrix of EQN. 27) into a vector and combining the vectors to make a single row to be used in the top level model. In this embodiment, step 1 transforms into the following vector:

(k_(1,1,1), k_(1,2,1) . . . k_(1,4,1), k_(1,1,2) . . . k_(1,4,2), k_(1,1,3) . . . k_(1,4,10))   EQN. 30

where the subscripts k_(x,y,z) are

x=process step

y=parameter

z=observation id

where “observation id” refers to a row of data from the individual process steps. Similarly, the other two steps would transform into the following two vectors

(k_(2,1,1), k_(2,2,1) . . . k_(2,4,1), k_(2,1,2) . . . k_(2,4,2), k_(2,1,3) . . . k_(2,4,10))   EQN. 31

and

(k_(3,1,1), k_(3,2,1) . . . k_(3,4,1), k_(3,1,2) . . . k_(3,4,2), k_(3,1,3) . . . k_(3,4,10))   EQN. 32

In this embodiment, these vectors (EQNS. 30, 31, and 32) are then combined to form a single row representing a single wafer that would have the following form:

(k_(1,1,1), k_(1,2,1) . . . k_(1,4,10), k_(2,1,1), k_(2,2,1) . . . k_(2,4,10), k_(3,1,1), k_(3,2,1) . . . k_(3,4,10))   EQN. 33

In this embodiment, the top level model would have the form of a matrix having as many rows as there are wafers. Each row would have the form of EQN. 33. Unfolded matrix (e.g., which can include a row represented by EQN. 33) can be used to solve EQN. 10 (e.g., the X portion of the PLS model). After X is solved, and Y is solved, the correlation between X and Y can be found (e.g., as described in the PLS and PSLDA method section above) to generate the top level model.

Method 2—Unfolding of Scores

Observation level data is unfolded into the top level model. One exemplary method uses raw data (e.g., as described above) and another method uses scores. In this embodiment, the MVA model is run on the lower level data as defined above in EQNS. 27-29. By way of example only, step 1 can have 2 principal components, step 2 can have 3 principal components, and step 3 can have 2 principal components. The scores matrices for these steps would have the following form, where the score matrices have as many columns as there are principal components for the step:

$\begin{matrix} {{{Step}\mspace{14mu} 1}:\begin{pmatrix} t_{1,1} & t_{2,1} \\ t_{1,2} & t_{2,2} \\ t_{1,3} & t_{2,3} \\ t_{1,4} & t_{2,4} \\ t_{1,5} & t_{2,5} \\ t_{1,6} & t_{2,6} \\ t_{1,7} & t_{2,7} \\ t_{1,8} & t_{2,8} \\ t_{1,9} & t_{2,9} \\ t_{1,10} & t_{2,10} \end{pmatrix}} & {{EQN}.\mspace{14mu} 34} \\ {{{Step}\mspace{14mu} 2}:\left( {\begin{matrix} t_{1,1} & t_{2,1} \\ t_{1,2} & t_{2,2} \\ t_{1,3} & t_{2,3} \\ t_{1,4} & t_{2,4} \\ t_{1,5} & t_{2,5} \\ t_{1,6} & t_{2,6} \\ t_{1,7} & t_{2,7} \\ t_{1,8} & t_{2,8} \\ t_{1,9} & t_{2,9} \\ t_{1,10} & t_{2,10} \end{matrix}\begin{matrix} t_{3,1} \\ t_{3,2} \\ t_{3,3} \\ t_{3,4} \\ t_{3,5} \\ t_{3,6} \\ t_{3,7} \\ t_{3,8} \\ t_{3,9} \\ t_{3,10} \end{matrix}} \right)} & {{EQN}.\mspace{14mu} 35} \\ {{{Step}\mspace{14mu} 3}:\begin{pmatrix} t_{1,1} & t_{2,1} \\ t_{1,2} & t_{2,2} \\ t_{1,3} & t_{2,3} \\ t_{1,4} & t_{2,4} \\ t_{1,5} & t_{2,5} \\ t_{1,6} & t_{2,6} \\ t_{1,7} & t_{2,7} \\ t_{1,8} & t_{2,8} \\ t_{1,9} & t_{2,9} \\ t_{1,10} & t_{2,10} \end{pmatrix}} & {{EQN}.\mspace{14mu} 36} \end{matrix}$

Alternatively, in another embodiment actual process data (e.g., see before) is used. A similar process as used above with the raw data method would translate the three matrices (e.g., EQNS 34-36) using scores, into the following vectors:

for step 1,

(t_(1,1), t_(2,1), t_(1,2), t_(2,2), . . . , t_(1,10), t_(2,10))   EQN. 37

for step 2,

(t_(1,1), t_(2,1), t_(3,1), t_(1,2), t_(2,2), . . . , t_(2,10), t_(3,10))   EQN. 38

and for step 3

(t_(1,1), t_(2,1), t_(1,2), t_(2,2), . . . , t_(1,10), t_(2,10))   EQN. 39

In this embodiment, the three steps would be combined into a single row of the top level model representing a single wafer and would look like

(t_(1,1,1), t_(1,2,1), t_(1,1,2), . . . t_(1,1,10), t_(1,2,10), t_(2,1,1), t_(2,2,1), t_(2,3,1) . . . t_(2,1,10), t_(2,2,10), t_(2,3,10), t_(3,1,1), t_(3,2,1) . . . , t_(3,1,10), t_(3,2,10))   (EQN. 40

where the subscripts x,y,z represent:

x=step

y=principal component

z=sample

Unfolded matrix (e.g., which can include a row represented by EQN. 40) can be used to solve equation 10 (the X portion of the PLS model). After X is solved, and Y is solved, the correlation between X and Y can be found (as described in the PLS section above) to generate a top level model.

Unfolding Multiple Processes

In some embodiments, unfolding multiple processes can work in a substantially similar fashion as unfolding multiple steps. In some embodiments, the method can include properly identifying wafers so the process data process to process remains within the same wafer.

In this embodiment, the processes can be defined as the steps are defined. In this embodiment, the lower level to top level unfolding has already been performed. The process matrices can look like:

$\begin{matrix} \begin{pmatrix} {t_{1,1,1,1}t_{1,1,2,1}t_{1,1,3,1}\mspace{14mu} \ldots \mspace{14mu} t_{1,s,v,n}} \\ {t_{2,1,1,1}t_{2,1,2,1}t_{2,1,3,1}\mspace{14mu} \ldots \mspace{14mu} t_{2,s,v,n}} \\ \vdots \\ {t_{w,1,1,1}t_{w,1,2,1}t_{w,1,3,1}\mspace{14mu} \ldots \mspace{14mu} t_{w,s,v,n}} \end{pmatrix} & {{EQN}.\mspace{14mu} 41} \end{matrix}$

where the subscripts w,s,v,n represent

w=wafer

s=process step

v=principal component (or alternatively, if modeling with raw data, variable)

n=observation number

In this embodiment, when multiple processes are combined, then the combined matrix will be structured according to:

$\begin{matrix} \begin{pmatrix} {t_{1,1,1,1,1}t_{1,1,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{1,1,s,v,n}t_{1,2,1,1,1}t_{1,2,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{1,2,s,v,n}\mspace{14mu} \ldots \mspace{14mu} t_{1,p,1,1,1}t_{1,p,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{1,p,s,v,n}} \\ {t_{2,1,1,1,1}t_{2,1,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{2,1,s,v,n}t_{2,2,1,1,1}t_{2,2,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{2,2,s,v,n}\mspace{14mu} \ldots \mspace{14mu} t_{2,p,1,1,1}t_{2,p,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{2,p,s,v,n}} \\ \vdots \\ {t_{w,1,1,1,1}t_{w,1,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{w,1,s,v,n}t_{w,2,1,1,1}t_{w,2,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{w,2,s,v,n}\mspace{14mu} \ldots \mspace{14mu} t_{w,p,1,1,1}t_{w,p,1,2,1}\mspace{14mu} \ldots \mspace{14mu} t_{w,p,s,v,n}} \end{pmatrix} & {{EQN}.\mspace{14mu} 42} \end{matrix}$

where the subscripts now account for

w=wafer

p=process

s=process step

v=principal component (or, alternatively, if modeling with raw data, variable)

n=observation number

EQNS. 27-42 as described above, define exemplary methods for combining multiple processes and steps. EQNS. 27-42 generate an “X Matrix” (e.g., an “nth” mathematical model) that can be used to solve EQN. 10 as described above. EQN. 10 can be combined with the “Y data” (e.g., metrology data) to complete a PLS top level model. In embodiments in which the models are based on scores rather than raw data, a top level matrix is generated that combines each of the multiple processes by using individual PLS models.

FIG. 2 is a schematic drawing depicting a method for implementing a model 200 for monitoring a manufacturing process for manufacturing a lot of semiconductor wafers 104, according to an illustrative embodiment of the invention. The method 200 involves performing at least two process steps 108A and 108B and acquiring data 204A and 204B (generally 204) for the output of the two process steps 108A and 108B, respectively. In this embodiment, the method 200 includes up to “n” process steps 108A, 108B, 108C and 108D. The data acquired (204A, 204B, 204C and 204D) for each of the “n” process steps is provided to a model (e.g., the model 124 of FIG. 1). The model 124 is then used to generate an output 128 that is a prediction of the occurrence of a fault and/or classification of the fault identified in one or more of the wafers. The model can be used to identify a fault or classify a fault in one or more of the process steps and/or in the final wafers output by performing each of the “n” process steps 108A, 108B, 108C and 108D.

FIG. 3 is a graph comparing a predicted yield 136 of a manufacturing process to an observed yield 140 of a manufacturing process, according to an illustrative embodiment of the invention. The predicted yield was created using a top level mathematical model, for example, the top level mathematical model 124 of FIGS. 1 and 2. The top level model was created using metrology data and process data from a first lot of 11 wafers (137 in FIG. 3). The model was then used to predict the yield for a second lot of 25 wafers (labeled as the prediction set 138). The graph in FIG. 3 displays two lines for the prediction set. The first line, labeled 136, is the predicted yield, based on the model. The second line, labeled 140, is the actual measured yield, measured after the completion of the device fabrication of the manufacturing process that produces the devices from the wafers (e.g., semiconductor devices). The similarity of the lines 136 and 140 demonstrates the efficacy of the model to predict process yield. The y-axis 148 represents yield of the process measured from 0, meaning no functional devices, to 1, meaning 100% functional devices on a given wafer.

The above-described techniques can be implemented in digital electronic circuitry, or in computer hardware, firmware, software, or in combinations of them. The implementation can be as a computer program product, e.g., a computer program tangibly embodied in an information carrier, e.g., in a machine-readable storage device or in a propagated signal, for execution by, or to control the operation of, data processing apparatus, e.g., a programmable processor, a computer, or multiple computers. A computer program can be written in any form of programming language, including compiled or interpreted languages, and it can be deployed in any form, including as a stand-alone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program can be deployed to be executed on one computer or on multiple computers at one site or distributed across multiple sites and interconnected by a communication network.

Method steps can be performed by one or more programmable processors executing a computer program to perform functions of the invention by operating on input data and generating output. Method steps can also be performed by, and apparatus can be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application-specific integrated circuit). Modules can refer to portions of the computer program and/or the processor/special circuitry that implements that functionality.

Processors suitable for the execution of a computer program include, by way of example, both general and special purpose microprocessors, and any one or more processors of any kind of digital computer. Generally, a processor receives instructions and data from a read-only memory or a random access memory or both. The essential elements of a computer are a processor for executing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto-optical disks, or optical disks. Data transmission and instructions can also occur over a communications network. Information carriers suitable for embodying computer program instructions and data include all forms of non-volatile memory, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto-optical disks; and CD-ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in special purpose logic circuitry.

The terms “module” and “function,” as used herein, mean, but are not limited to, a software or hardware component which performs certain tasks. A module may advantageously be configured to reside on addressable storage medium and configured to execute on one or more processors. A module may be fully or partially implemented with a general purpose integrated circuit (“IC”), FPGA, or ASIC. Thus, a module may include, by way of example, components, such as software components, object-oriented software components, class components and task components, processes, functions, attributes, procedures, subroutines, segments of program code, drivers, firmware, microcode, circuitry, data, databases, data structures, tables, arrays, and variables. The functionality provided for in the components and modules may be combined into fewer components and modules or further separated into additional components and modules. To provide for interaction with a user, the above described techniques can be implemented on a computer having a display device, e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor, for displaying information to the user and a keyboard and a pointing device, e.g., a mouse or a trackball, by which the user can provide input to the computer (e.g., interact with a user interface element). Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback, e.g., visual feedback, auditory feedback, or tactile feedback; and input from the user can be received in any form, including acoustic, speech, or tactile input.

The above described techniques can be implemented in a distributed computing system that includes a back-end component, e.g., as a data server, and/or a middleware component, e.g., an application server, and/or a front-end component, e.g., a client computer having a graphical user interface and/or a Web browser through which a user can interact with an example implementation, or any combination of such back-end, middleware, or front-end components. The components of the system can be interconnected by any form or medium of digital data communications, e.g., a communications network. Examples of communications networks, also referred to as communications channels, include a local area network (“LAN”) and a wide area network (“WAN”), e.g., the Internet, and include both wired and wireless networks. In some examples, communications networks can feature virtual networks or sub-networks such as a virtual local area network (“VLAN”). Unless clearly indicated otherwise, communications networks can also include all or a portion of the publicly-switched telephone network (“PSTN”), for example, a portion owned by a specific carrier.

While the invention has been particularly shown and described with reference to specific illustrative embodiments, it should be understood that various changes in form and detail may be made without departing from the spirit and scope of the invention. 

1. A system for creating a mathematical model according to a multivariate analysis, the system comprising: a. means for acquiring metrology data for semiconductor wafers at the conclusion of a final process step for a manufacturing process; b. means for acquiring data for a plurality of process variables for a first process step for manufacturing semiconductor wafers; c. a model generation module for creating a first mathematical model of the first process step based on the metrology data and the acquired data for the plurality of process variables for the first process step; d. means for acquiring data for a plurality of process variables for at least a second process step for manufacturing the semiconductor wafers and a model generation module for creating a second mathematical model of the first process step based on the metrology data and the acquired data for the plurality of process variables for the first process step; e. a model generation module for creating an n^(th) mathematical model based on the metrology data and the data for the plurality of process variables for each of the n process steps; and f. a model generation module for creating a top level mathematical model based on the metrology data and the models created by elements c., d. and e., the top level mathematical model being based on those process variables that have a substantial effect on the metrology data.
 2. The system of claim 1, comprising a processor for comparing data for each process step for manufacturing a semiconductor wafer to the top level model generated by element f. to determine whether a fault occurred during any of the n process steps.
 3. The system of claim 1, comprising a processor output connected to at least one tool used to perform at least one process step.
 4. The system of claim 3, wherein the processor modifies an operating property of at least one process step based on whether the top level model determines that a fault has occurred during at least one process step.
 5. The system of claim 3, wherein the processor identifies which of the process variables causes manufacturing yield to be below a predetermined limit.
 6. The system of claim 1, comprising a processor for comparing data for the first process step for manufacturing a semiconductor wafer to the top level model generated by element f. to determine whether manufacturing yield will be above a predetermined limit.
 7. The system of claim 1, wherein the process variables that have a substantial effect on the metrology data are weighted more in creating the top level mathematical model than process variables that do not have a substantial effect on the metrology data.
 8. The system of claim 1, wherein the top level mathematical model is a yield prediction model based on the process variables selected by the models created by elements c., d. and e. 